d=read.table("clipboard",header = T)
source("F:/msaR.R")
D=msa.X(d)
Y=scale(D,center=F,scale=T) #标准化
View(Y)
#1确认数据是否是否适合做因子分析 有kmo和bartlett进行验证
msa.KMO(Y)
msa.bartlett(Y)
#2构造因子变量
FO=msa.fa(Y,3,rotation = "none")#FO未旋转分为三个因素
FO
#3旋转因子更具可解释性
FR=msa.fa(Y,3,rotation = "varimax")#旋转
FR
#4计算因子得分并作因子图
biplot(FR$scores[,1:2],FR$loadings[,1:2])
biplot(FR$scores[,c(1,3)],FR$loadings[,c(1,3)])
biplot(FR$scores[,2:3],FR$loadings[,c(2,3)])

